Thermally excited spin-current in metals with embedded ferromagnetic nanoclusters 
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We show that a thermaUy excited spin-current naturally appears in metals with embedded ferro- 
magnetic nanoclusters. When such materials are subjected to a magnetic field, a spin current can 
be generated by a temperature gradient across the sample as a signature of electron-hole symmetry 
breaking in a metal due to the electron spin-flip scattering from polarised magnetic moments. Such 
a spin current can be observed via a giant magneto-thermopower which tracks the polarisation state 
of the magnetic subsystem and is proportional to the magnetoresistance. Our theory explains the 
recent experiment on Co clusters in copper by S. Serrano-Guisan et al [Nature Materials AOP, 
doi:10.1038/nmatl713 (2006)] 
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Theoretical studies of spin-currents in metals and semi- 
conductors have recently been fuelled by the search for 
materials to develop spintronic devices [l|. Among the 
materials studied of particular interest had been com- 
posite ferromagnetic and nonmagnetic metallic systems, 
such as metallic multilayers 0] and metals with embed- 
ded ferromagnetic nanoclusters where the electronic 
transport strongly depends on the electron spin leading 
to the effect of the giant magnetoresistance (MR). 

In metals with ferromagnetic nanoclusters 0i B or 
magnetic impurities j5|] with a large spin, s ^ 1, the mag- 
netoresistance reflects the degree of polarisation, {1^) = 
coth (y) — y~^, y = fiBs/kT of a magnetic sub-system by 
a magnetic field B = n^B, 



i?(0) 
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The MR effect develops across the magnetic field range 
Bmb. ~ kT/fj,s, and its strength depends on the rela- 
tive difference between the mean free path time t^{i) of 
electrons with spins parallel(antiparallel) to the cluster 
polarisation axis. Here, T|(|) is defined as the mean free 
path in a metal with magnetic clusters whose spin polar- 
isations are fully up (down). The Fermi density of states 
per spin, 7, is identical for both electron spin orienta- 
tions, and r = ^ {t^ + t^) is a characteristic mean free 
path time, while /i is the Bohr magneton. 

In this Letter we show that a thermally excited spin- 
current naturally appears in metals with embedded fer- 
romagnetic nanoclusters (FmnC's) when subjected to a 
magnetic field and temperature gradient. This effect ap- 
pears through the manifestation of electron-hole symme- 
try breaking via electron spin-flip scattering from po- 
larised magnetic moments. When two parts of a metal- 
lic sample with embedded FmnC's are held at different 
temperatures Ti > T2, the thermally equilibrating heat 
flux must be accompanied by a transport of magnetisa- 
tion 00,13) a spin current, jg. This spin current 
acts to equilibrate the cluster polarisation on opposing 
sides of the temperature gradient. Locally, the equilibra- 
tion process requires spin transfer from clusters to con- 
duction electrons via spin- flip scattering at a rate t^^, 



which is the inverse mean free path time of the spin-flip 
scattering process (formally defined latter, Eq.lO). After 
a single spin- flip scattering event, the scattered electron 
will carry the transferred magnetisation while diffusing 
(with diffusion coefficient D) between different parts of 
the sample held at different temperatures, thus leading 
to a spin-current. 



=/3Vr, /3; 



(1) 



Here / = 



and X = fiB/kT show that the spin 



current will persist up to a value of the magnetic field 
Bg ~ kT/ 11, which is much larger than the typical field 
value at which MR develops, Bmr <C Bs- Additionally, 
due to the difference in the mean free path times of "up" 
and "down" spin carriers, (r^ — t^) (Iz), the spin current 
'drags' a charge current, j cx xVT . In an open circuit, 
this generates a thermopower, V12 cx c {B) [Ti ~T2], with 
a strong magnetic field dependence resembling (for B < 
kT/ii) that of MR, 



S(i?) ^ c(B) - c(0) 



e Ts (r^ - r|) 



A(i?). (2) 



A microscopic justification of the above-presented phe- 
nomenological argument is supported by the following 
analysis of the electron-hole asymmetry in a compos- 
ite metal with ferromagnetic components. In normal 
metals, the electron- hole asymmetry leading to the ther- 
mopower is caused by the energy dependence of the den- 
sity of states of electrons near the Fermi level. In ma- 
terials containing a polarised magnetic subsystem, the 
electron-hole asymmetry is created in an alternative way 
— via the enerry and spin dependence of a quasi-particle 
scattering rate [a, 0, ■ In the system we discuss in this 
Letter, the formation of kinetic electron-hole asymmetry 
can be illustrated using Kondo-type model which 
treats a FmnC as an impurity with a large spin, s, and 
contains all the necessary ingredients to describe the MR 
effect simultaneously with the generation of spin-current 
and magneto-thermopower. 



2 



Thus, we model the FmnC's by the Hamiltonian 



a) 



b) 



((S-')2)G-' 



Here, the potential U accounts for the FnmC charge and 
for band mismatch between normal and magnetic metals, 
whereas J is the exchange interaction; /is is the magnetic 
moment of each cluster; a stands for the electron spin 
operator. The MR analysis in materials with large spin 
clusters, s 3> 1 can be done using a static exchange field 
model |3|, where the operator s is replaced by si (here, 1 
is a unit vector in the direction of polarisation of an in- 
dividual cluster). In contrast, the analysis of a thermally 
excited spin current and the magneto-thermopower needs 
to take into account the quantum nature of the cluster 
spin, manifested through the electron spin-flip process. 

The formation of electron-hole asymmetry in the quasi- 
particle lifetime becomes apparent after analyzing the 
imaginary part of the self-energy to the lowest non- 
vanishing order of perturbation theory. Fig. ^ shows 
the Keldysh diagrams which appear in the self-energy 
j^A{R) ^YiQ second order of the electron-cluster in- 
teraction, after averaging over random positions, r^, of 
the clusters. We assume a thermal distribution of clus- 
ter spin polarisation, with temperature T. The dia- 
grams in Fig.^incorporate two types of vertices: crosses 
stand for the potential U, and dots for the exchange 
coupling J. Dashed lines represent impurity averaging 
[ = ncS{r — r')]. Thin solid lines indicate the thermo- 
dynamic average and dynamical correlators of the cluster 
spin: a loop stands for the thermodynamic average (sz), 
a double line without an arrow is ((s^)^), whereas a thin 
double line with arrows stands for the free-spin Green 
bmctions 10], D^/^{t) = ^i{[s+{t)s- {0)]^)e{±t) and 

= li({s+(t)s-(0)}^), where s± 
the energy representation. 



Sx ± isy. In 



uo — fiB ± iS 

i?^-M = coth^(z?«H-2?^H). (3) 

Although Wick's theorem cannot be applied to spin op- 
erators in higher orders of perturbation theory, this does 
not affect the analysis of the diagrams in the Born ap- 
proximation. Finally, bold solid lines with arrows are free 
electronic Green functions. 
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e — £,p ± iS' 



G^{e, p) = tanh — (G«(e, p) - G^(e, p)) , (4) 

where the electron energy e is determined with respect 
to the Fermi level. 

The diagram in Fig. ^ (a) describes scattering from a 
scalar potential, whereas diagrams (b) and (d) only in- 
clude the exchange interaction. The interference between 
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JS-a+-''jS+a- qAijK JS+a-^''jS-a+ ^0^^ 

FIG. 1: Self-energy digrams: thick solid line is free electrons 
Green's function, G± — G{e ± /i-B), thin dashed line is im- 
purity line, thin solid line connected both ends to a single 
vertex is an average of spin operator, double thin solid line 
connected to two vertices is the spin-spin correlator, cross is 
the potential vertex U and dot is the exchange vertex J. 



the scalar and exchange potential scattering amplitudes 
(from the same impurity) is taken into account by the 
diagrams in Fig. 1(c), which is responsible for the spin- 
dependence of the electron mean free path in the pres- 
ence of polarised clusters. Together, these diagrams yield 

mm 

ImS^(£) = 7r7nc [U^ + J^s{s + 1) + 2UJ{s,^)a,] - 

- 7r7nc J^(s3)aztanh (^-^-^f-^^ > (5) 

where ric is the concentration of FmnC's. Diagrams (a)- 

(c) describe elastic processes. The spin-flip diagrams 

(d) contain an inelastic part resulting in the energy- 
dependent contribution towards ImS^ in Eq.®, and 
they are responsible for the kinetic electron-hole asym- 
metry. This asymmetry is most pronounced when the 
ensemble of clusters is fully polarised, (Iz) 1, and van- 
ishes when (Iz) = 0. 

To explain the origin of this asymmetry and its polar- 
isation dependence we consider the limit of (Iz) — s- 1, so 
that all FmnC's are polarised 'up'. Fig. 2(a) shows that 
an incident spin-| electron with energy e > fj,B above 
the Fermi level is able to flip its spin turning the spin of 
FmnC from the polarisation axis. The amplitude of such 
a process is A J{s — l|s^|s) ~ J\/s, which results in 
the scattering rate c>c J^s. The relevant range of a 
quasiparticle excitation energy for this process is set by 
the FmnC energy splitting, /iB, between its initial and 
scattered spin state, which differ by Sz — = 1. A sim- 
ilar process is possible for the incident spin-t hole with 
e < —fiB below the Fermi level. Fig. 2(d): it corresponds 
to the reverse process of an equilibrium spin-| electron 
relaxing into an empty state with opposite spin below the 
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holes \iB electrons 

FIG. 2: Schematic representation of allowed and forbidden 
processes leading to the energy and spin dependent quasi- 
particles relaxation rate ( a), d) - allowed processes, b), c) - 
forbidden processes, e) - dependence of Im E'* as function of 
energy for different spin directions - graph is drawn for the 
case of JU < 0) 



Fermi level (hole). In contrast, neither spin-^ electrons 
above the Fermi level, nor spin-| holes shown in Figs. 
2(b,c) can scatter with changing spin state, since they 
cannot increase further the maximal = s of a clus- 
ter. All together, these processes determine the energy- 
dependent part in the electron scattering rate from fully 
polarized clusters, 



ImE 
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'^T(i) ^ 



T„ = ; J S. 



([/^ + J^s(s + 1) ± 2[/Js), 



(6) 



Note that for a cluster with N magnetic atoms S oc N 
and U (X N so that r^^J^ cx N^, whereas cx N. The 
resulting behaviour of the quasi-particle lifetime for elec- 
trons in the presence of polarised FmnC's is sketched in 
Fig. Ele). It indicates that the electron- hole asymmetry 
is inverted in the opposite spin channels, in accordance 
with the symmetry of the Kondo problem [l^. 

The inverted electron-hole asymmetry in the opposite 
spin channel results in the generation of a spin current 
when two parts of the system are held at different temper- 
atures. To describe a thermally generated spin-current, 
js = /3VT, as well as the electric current j = xVT, we 
use the kinetic equation approach. That is, we study a 



steady-state kinetic equation. 



[v • V - eE • ^p] g± = {I±{g, s^)) , 



(7) 



where the collision term I±{g,Sz) describes the balance 
in the distribution functions of spin-f and p_ of spin-| 
electrons which scatter from a group of FmnC's in a given 
initial spin state, Sz- The brackets (...) stand for aver- 
aging over the thermal distribution of cluster spins. The 
averaged collision term therefore takes into account all 
electron scattering processes described by the diagrams 
in Fig. HI 



2TTnc\U±Jsz 



J {q'± ~ 9±) S i^p - ep') 
si + s±Sz) g'^ (1 - g±) 



-{s^^sl+sT Sz) g± (1 - g'^)] S {ep, - Ep T fJ-B) . 

The first term in I± describes the elastic spin-conserving 
processes in Fig. ffl^a-c), whereas the second takes 
care of spin-flip processes corresponding to diagrams in 
Fig. ^d). By inspection, one can see that the probabil- 
ities of electrons i ^ t and t ^ i spin-flip scattering 
from the same cluster with a given Sz slightly differ. In 
a system where FmnC's are partially polarised (along 
the external magnetic field B = n^i?) this leads to the 
energy- and spin-dependent electron momentum relax- 
ation rate, 2lm'E'^^/h (where a = t,|) given by Eq. 

. In the leading order in a small temperature gradient 
(vftT'^VT < 1), the spin-current, js = | / 7deTr{(T^p} 
and electric current, j =f / 7deTr{p} are determined by 
the first angular harmonic in momentum space, p of the 
electron spin-density matrix, ^ « gx + p ■ p/p- This first 
angular harmonic can be found using the equation 



^ImS^p: 

h 



evpd^gT 
f 



VT, 



where gris) = [e'^/'^^ + l] "'^ is the Fermi function. 

As a result, we calculate the spin-current, js = /3VT, 
which, for large spin clusters, s 3> 1 is described by 



f3- 



(8) 



(n + n) (U' 



ir, + nf-{T,-nf{iz)^ 



2 • 



and 



with 



In the above result, scattering rates ''"ij-) 
defined in Eq. ©, and the factor / = i^^^^iyi 
X = fiB/kT is specified for the case s ^ 1 and takes 
into account the suppression of the spin transfer rate 
at low temperatures and high magnetic field, such that 
fiB > kT, since it requires the electron energy trans- 
fer e — e' = iiB. Also, in most of the metals with em- 
bedded ferromagnetic clusters, the maximum MR effect 
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is A < 10%, which is additionally suppressed by non- 
magnetic impurities and phonon scattering. Therefore, 
we can simplify the result in Eq. ||SJ) further, using the 
fact that (t| — T| ) / (t| + T| ) <C 1 , which leads to the 
approximate form in Eq. 

Due to the difference, (r| — r|) (Z^) between the mean 
free path times of spin-l and -J, carriers scattering from 
the ensemble of partially polarised FmnC's, the spin- 
current 'drags' a charge current j (x >c\7T. For s 1, we 
find 



3t, 



(n + Tif - (n -nf i^zf 



Together with the MR of this material obtained using the 
same approximations. 



R 



3 n + n 



1 



(n 



(n + 



ih) 



this determines the thermopower coefficient c = —kR 
with a magnetic field dependence. The magneto- 
thermopower = c{B) — c(0) is related to the ob- 

servable MR as 



A(B) 



er,(rx-r|) 1 + A(B) 



/ 



(9) 



The magnetic field dependence of implicit in the 

above result, contains two field scales [15j- At a low 
magnetic field where the polarisation of clusters devel- 
ops, the MTP is proportional to the MR and saturates 
together with A(i3) at the field i?MR- At a higher 
field range, B > kT/fi the MTP is suppressed and 
dies away. For weak MR materials, A < 0.1, where 
T| — Tf ^ T = ^ {ti + T^) , the above relation can be 
approximated using the formula in Eq. (0). 

At the early stage of study of Kondo impurity in 
metals ^3 a strong magneto-thermopower has been no- 
ticed in various dilute magnetic alloys 0, 0, 0, ^| . 
Later, a strong magnetic field dependent thermopower 
has been observed in metals with embedded ferromag- 
netic nanoclusters USUI mill 113, where it has been 

attributed to the energy dependent density of states in a 
weakly ferromagnetic metal with a complex band struc- 
ture [9j. However, when viewed in terms of the above- 
presented theory, these observations indicate the pres- 
ence of thermally excited spin-current, Eqs. (|1I8I) gener- 
ated by a temperature gradient. 

Recently, a giant MTP has been measured in copper 
with embedded Co clusters of a controlled size. It was 
noticed that MTP weakens with increase of the num- 
ber of atoms, TV, in a FmnCf2^. This behaviour as 
well as the observed direct correspondence between MTP 
and MR can be explained on the basis of the result in 
Eq.Q. The matter is that for a cluster with N mag- 
netic atoms T^~^J^ (X 7V^, whereas (x N. Therefore, 



for two samples with the same densities of the embed- 
ded clusters containing N and N' atoms, we estimate 
E{N)/E{N') cx N'/N, the MTP is weaker in the system 
with larger clusters. 
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